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Electrochemical setup — a unique chance to simultaneously control 
orbital energies and vibrational properties of single-molecule junc¬ 
tions with unprecedented efficiency! 
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Abstract: 

Impressive advances in nanoscience permit nowadays to manipulate single molecules and broadly control many of their prop¬ 
erties. Still, tuning the molecular charge and vibrational properties of single molecules embedded in nanojunctions in broad 
ranges escaped so far to an efficient control. By combining theoretical results with recent experimental data, we show that, 
under electrochemical control, it is possible to continuously drive a redox molecule (viologen) between almost perfect oxidized 
and reduced states. This yields an unprecedentedly efficient control on both vibrational frequencies and the surface-enhanced 
Raman scattering (SERS) intensities. The broad tuning achieved under electrochemical control by varying the overpotential 
(“gate potential”) within experimentally accessible ranges contrasts to the case of two-terminal setups that require high biases, 
which real nanojunctions cannot withstand. The present study aim at stimulating concurrent transport and SERS measurements 
in electrochemical setup. This may open a new avenue of research that is not accessible via two-terminal approaches for better 
understanding the transport at nanoscale. 


Keywords: nanotransport; single-molecule junctions; electrochemical scanning tunneling microscopy; surface-enhanced Raman 
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1 Introduction on single-molecule junctions represent significant attempts in 

this direction. Since the first demonstration that SERS can 


Despite important advances in the last decades, molecular 
electronics remains confronted with a series of difficulties. 
Many of them result from unsatisfactory characterization of 
a molecule under in situ conditions. Properties of a molecule 
embedded in a nanojunction in a current-carrying state may 
differ from its equilibrium properties in a manner reminding 
differences between biological cells in in vivo and in vitro 
situations. Eundamental processes at nanoscale and tailoring 
molecular devices for practical applications can be better un¬ 
derstood if transport data can be correlated with other molec¬ 
ular properties obtained from independent measurements of 
a different kind. Vibrational properties deduced via surface- 
enhanced Raman spectroscopy (SERS) belong to this cate¬ 
gory, since they can provide valuable complementary infor¬ 
mation needed for a better in situ characterization. 

Recent concurrent SERS— and ohmic conductance studies 
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be employed to study molecular junctions,— there is a con¬ 
tinuing interest to utilize this technique in molecular elec¬ 
tronics.—"— The electrodes of molecular junctions can act as 
highly efficient plasmonic antennas.— The enormous local 
electric field (especially close to the sharp tip of a scanning 
tunneling microscope-i or a nanoparticle^) can have a dra¬ 
matic effect on the surface-enhanced Raman scattering from 
molecules in junctions. Correlating SERS data with simul¬ 
taneously acquired conductance data provides important evi¬ 
dence on the chemical identity of the active molecule, on how 
it bonds to electrodes, and on experimental conditions {e.g., 
solvent, sample treatment, etc). At low biases corresponding 
to a linear transport regime, these refer to properties of the 
molecule (linked to electrodes but) at equilibrium (fluctuation- 
dissipation theorem). In particular, they refer to a molecule in 
a given charge (often neutral) state. 

Applying higher biases I4 on ^ molecular junction may 
change molecular properties, also including the charge of the 
molecule, and this change can reflect itself in a change of 
its vibrational properties. A recent SERS study on fullerene- 
based electromigrated junctions^ found that vibrational fre¬ 
quencies are significantly shifted under applied bias. Com- 
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panion density functional theory (DFT) calculations indicated 
that the observed frequency shifts are inconsistent with a sim¬ 
ple vibrational Stark effect, but they can result from a bias- 
driven change of the electronic charge of the molecule.— By 
applying realistic source-drain voltages I 4 in a two-terminal 
setup, it is only a very partial reduction (in contradistinction 
to a complete reduction of the molecular species, correspond¬ 
ing to changing the molecular charge by an entire electron) 
that can be achieved (c/ Ref. [12 and Fig. Si. 

The theoretical results reported in this paper, obtained 
by combining a model study backed by companion quan¬ 
tum chemical calculations and existing experimental data,— 
demonstrate that an SERS study in electrochemical scan¬ 
ning tunneling microscopy (EC-STM) setup can provide valu¬ 
able information complementary to that obtained via a single¬ 
molecule transport study. The EC-STM approach exploits the 
flexibility of a three-terminal setup: both the bias I 4 = V) —14 
between the STM-tip (f) and substrate (s) and the overpoten¬ 
tial 77 = Veq — can be independently controlled. Eor 

viologen-based junctions,— the equilibrium potential Vgq — 
—0.46 V. Charge transport through several redox systems in 
EC-STM setup has been experimental studie d but cor¬ 

responding SERS experiments are missing. The advantages 
of this setup introduced in Tao’s seminal work,— which en¬ 
ables a practically complete reduction/oxidation of the molec¬ 
ular species, for concurrent SERS and transport studies will 
be emphasized. 

2 Methods 

To demonstrate the usefulness of a concurrent SERS and trans¬ 
port study through a redox unit embedded in a single-molecule 
junction in a three-terminal electrochemical setup, we will 
consider viologen-based molecular junctions, which were in¬ 
vestigated in a comprehensive experimental study.— 

The core of this molecule consists of a redox-active 4,4'- 
bipyridinum dication (44BPY++). The parent (neutral bipyri¬ 
dine 44BPY*’) molecule embodied in nanojunctions formed 
the object of numerous transport, transport-related—"— 
and SERS studies. This redox-active molecule was 
used, e.g., as backbone in self-assembled monolayers^— and 
in various functional materials .2ii^ It is a showcase redox 
molecule; the first oxidation-reduction process 44BPY++ ^ 
44BPY+*, which will be examined below, is completely re¬ 
versible in bulk solutions.— 

High-level quantum chemical calculations at the den¬ 
sity functional theory (DET) level using the Becke’s three- 
parameter hybrid functional B3LYP and basis sets of triple- 
zeta quality augmented with diffuse functions (aug-cc-pVTZ) 
as implemented in the GAUSSIAN 09 package^i^ have been 
performed for geometry optimizations and for obtaining the 
Raman spectra of the oxidized and reduced viologen core (di¬ 


cation 44BPY++ and radical cation 44BPY+*, respectively), 
which are the charge species that contribute to the measured 
current— in the specific case considered in this paper. 

In addition to the validation against the experimental trans¬ 
port data, we have also microscopically validate the Newns- 
Anderson model employed below by calculating the low¬ 
est electron affinities at the EOM-CCSD (equation-of-motion 
coupled cluster singles and doubles) level,— which represents 
the quantum chemistry state-of-the-art for molecules of this 
size; see the ESIf. 

To compute the vibrational frequencies and Raman scatter¬ 
ing intensities of a molecular junction under bias (V), ^ 0 ,77 7 ^ 
0 ), we have resorted to an interpolation weighting method de¬ 
scribed in the ESIj, wherein the weight is expressed in terms 
of the Vb- and 77 -dependent LUMO occupancy. 


3 Results 

The Raman spectra of the dication (44BPY++) and cation 
(44BPY+*) species in acetonitrile computed as described in 
ESIf are presented in Table S1 and Fig. [T] Fig.El Fig. SI, 
Fig. S2, Fig. S3, and Fig. S4. (Throughout, label S refers to 
ESIf). These results reveal notable differences between vi¬ 
brational properties of the two different redox charge species. 
Both vibrational frequencies and Raman scattering intensities 
significantly depend on the charging state. 
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Fig. 1 Raman scattering activities (spectral lines in A'^/a.m.u.) of 
the dication 44BPY++ in acetonitrile. The envelope (red line) has 
been obtained by convoluting the computed spectral lines (green 
spikes) with Lorentzian functions of half-width 20 cm"*. 

For SERS observability, the important issue is whether 
a transport setup permits to broadly control the molecular 
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Fig. 2 Raman scattering activities (spectral lines in A^/a.m.u.) of 
the cation 44BPY+* in acetonitrile. The envelope (red line) has been 
obtained by convoluting the computed spectral lines (green spikes) 
with Lorentzian functions of half-width 20 cm^*. 


charge, enabling a continuously switching between the dica- 
tionic (44BPY++) and cationic (44BPY+*) species. They cor¬ 
respond to an oxidized (n/ = 0) and a reduced («/ = 1) LUMO, 
respectively. As a central point of the present analysis, we 
have used the transport data of Ref. [ij to show that the EC- 
STM transport setup does enable this switching. 

In an EC-STM setup, transport data can be acquired in 
two basic modes; constant bias and variable bias modes. In 
constant bias mode, the STM-tip (f) and substrate (s) poten¬ 
tials Vt,s are varied such that Vh =Vt —14 is kept constant. 
In variable mode, I 4 is varied at constant substrate potential 
14 (J? = const). As visible in Fig. [3 and Fig. @1 the theoreti¬ 
cal curves successfully reproduce the experimental currents^ 
measured in both aforementioned modes. Because the Newns- 
Anderson model utilized to obtain these theoretical curves has 
been discussed in detail elsewhere, only a few relevant details 
are given below and in the ESIf. 

An aspect worth to mention is the LUMO position. The 
LUMO energy utilized in the transport calculations lies at 
—eVeq — 0.46 eV above electrodes’ equilibrium Eermi energy, 
in agreement with the experimental data.— Noteworthy, this 
value, which implicitly enters the definition of the overpoten¬ 
tial (see Sec. [T] and, for example, the Supporting Information 
of ref. [ 3 ^ . corresponds to a LUMO energy of the embed¬ 
ded molecule, is different from that of the isolated molecule. 
Schemes to disentangle this energy difference in contributions 
with clear physical origin have been discussed recently. 

As another particularly relevant detail, we mention the sub¬ 
stantial asymmetry of the molecule-electrode couplings 5 = 


Ef /r, where E = (E^ + E,) /2, of the viologen-based EC-STM 
junctions. *^ The value 5 <C 1 (see legends of Fig.|4]and Fig.O, 
or alternatively E, < E„ realistically accounts for the experi¬ 
mental setup of ref. [O]; the viologen molecule is chemically 
bound to the substrate (s) but not to the STM tip (f). Albeit 
substantial, this asymmetry is not so pronounced as recently 
found for azurin-based EC-STM junctions.— 


Once validating the transport model, we have employed it 
to compute the bias dependent LUMO occupancy «/, which is 
also shown in Fig. |3] and Fig.|4] 



Fig. 3 The present theoretical model, described in detail 
elsewhere, ^^ 1 ^^ is able to reproduce the currents I measured by 
varying the overpotential 7] in constant bias mode for 
single-molecule junctions based on viologe n and to show that 
the LUMO occupancy ra; can be continuously tuned between 0 
and n; ~ 1, which correspond to almost perfect oxidized and 
reduced states, respectively. Relevant details on the model and the 
parameters <^, A, y, E, and S are given in ESIf. The experimental 
current (red curve) presented here was obtained by digitizing 
Fig. 8A of ref. [3, where the value of the preset current is 
ho = 0.1 nA; notice that Ej, = -f 0.05 V of ref. [T^ corresponds in the 
present notation to Vh = —0.05 V. (Currents scaled such that the 
maximum theoretical current is equal to unity.) 


As alternative to other approaches to SERS in biased molec¬ 
ular junctions,— to estimate the vibrational frequencies (Oy 
and the Raman scattering intensities Ay of the various modes 
V we adopt here an interpolation method described in ESIf. 
The Vh- and t]-dependencies of cOy and Ay follow from those 
of «/. They are depicted in Fig. 0 and Fig.|6l The dependence 
Hi = ni{ri) shown in Fig.[3]represents a key point of the present 
analysis. 
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Fig. 4 The Newns-Anderson model with reorganization, described 
in detail elsewhere, is able to reproduce the experimental curve 
1 — V curve measured for single-molecule junctions based on 
viologen.— In the bias ranged sampled in experiment— (depicted 
by the red experimental curve and indicated by the two black 
vertical dashed lines) the LUMO occupancy is negligible (n; ~ 0). 
This demonstrates that an effective reduction cannot be achieved by 
using source-drain voltages in the experimentally accessed V),-range. 
Relevant details on the model and the parameters y, F, and 8 
are given in ESIf. 


4 Discussion 

Fig. [3 and Fig. |4] show that, in principle, the charge of a 
molecule embedded in a biased EC-STM junction can be con¬ 
trolled both in constant bias mode and variable bias mode. 
However, as far as the bias ranges (that can be) accessed in 
experiment— are concerned, there is an important quantitative 
difference between the two operating modes. 

Within the whole V),-range that has been sampled in ex¬ 
periment (|14| < 0.4 V—) the LUMO occupancy is negligible 
(«/< 0 . 02 , c/ Fig.©. This corresponds to an almost perfect 
oxidized 44BPY++ state. As shown in Fig. |5l this yields a 
negligible 14 -dependence of the vibrational frequencies and 
Raman scattering intensities. A significant change in LUMO 
occupancy, also accompanied by significant variations of the 
vibrational properties {cf. Fig.|5]l, can only be achieved at sub¬ 
stantially larger 14 ’s, which molecular junctions could hardly 
withstand. This is a general feature of the off-resonant tunnel¬ 
ing in two-terminal setup. 

Fig. [3 depicts a totally different situation. By varying the 
overpotential 77 within the range accessed in experiment,— 
the LUMO occupancy can be continuously tuned between an 
almost perfect oxidized state 44BPY++) and an al- 



Fig. 5 Dependence on the source-drain bias I 4 of the changes in 
frequencies and Raman scattering activities (in A^/a.m.u.) of several 
representative vibrational modes specified in the legend. As visible, 
in the range sampled in experiment,— which is indicated by the two 
black vertical dashed lines, the vibrational properties do not 
significantly vary with I 4 . 

most perfect reduced state (u/A-l, 44BPY+*). Because con¬ 
stant bias experiments can practically sample the whole range 
0 < «/ < 1 , the values (Oy and Ay can continuously cover the 
values of the various (v) vibrational modes corresponding to 
the dicationic and cationic species. The complete list of cUy’s 
and Ay’s for 44BPY++ and 44BPY+* is presented in Table 
SI 77 - and 14 -dependencies of several representative vibra¬ 
tional modes are shown in Fig. |5] and Fig.| 6 ] In contrast to the 
insignificant impact of I 4 (Fig. ID, 77 -driven variations cUy = 
COy(ri) and Ay =Av{ri) like those depicted there are substan¬ 
tial (Fig. ID. Frequencies and intensities of the various Raman 
active modes are affected in different ways. For some modes, 
the state of charge only has a weak impact on the frequency 
while the intensity is strongly affected and vice versa. The fre- 
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Fig. 6 Dependence on the overpotential T] of the changes in 
frequencies and Raman scattering activities (in A^/a.m.u.) of several 
representative vibrational modes specified in the legend. 


quency of mode 10 (ring out-of-plane deformation) exhibits 
the highest frequency change; CUio.c — ®io,d = 157.6cm^'. 
Although the Raman activity decreases by one order of mag¬ 
nitude upon reduction {cf. Table SI), the small intensity of 
this mode would probably be challenging for experimental¬ 
ists. At the other extreme, the LUMO reduction only yields 
a decrease in frequency of mode 49 (CH stretch) amounting 
to ©49 c — ®49 D = 19.7 cm^'. However, the Raman intensity 
of this mode (which is the highest for the dication, cf. Table 
SI) is diminished by a factor A 49 £)/A 49 c « 22. This makes it 
a good candidate to be monitored in experiments. The other 
modes shown in Fig.|5]and Fig.| 6 ]are those identified in exper¬ 
imental Raman spectra of the radical anion 44BPY^*;— see 
ESIf for further details. Here we only note that mode 64 is 
related to the so-called quinoidal distortion, corresponding to 
a shortening of the inter-ring C-C bond and of the C-C bond 


parallel to it, and a lengthening of the C-C bond between them 
as well as of the C-N bond.^^ 


In an EC-STM setup similar to that of Ref. [131, the electro¬ 
magnetic (EM) enhancement occurs in a tiny region near the 
STM tip, and the charge transport mainly proceeds through 
the single (physisorbed) molecule which is closest to it. Still, 
other molecules nearby may also feel a significant electro¬ 
magnetic enhancement and may also experience LUMO en¬ 
ergy shifts driven by more or less similar local I 4 and 77 val¬ 
ues {cf. Eq. (SI)]; the spatial potential profile is hard to con¬ 
trol/determine. So, the recorded SERS dependence on 77 (and 
Vb) may not be (entirely) due to the (single) molecule mediat¬ 
ing the charge transport. Therefore, correlating the measured 
SERS with the currents measured by varying 77 (or Vb) rather 
than with 77 (or Vb) is preferable; it would be the most clear 
indication that the Raman signals come from the particular 
molecules inside the nanogap that are responsible for the cur¬ 
rent between the electrodes. Results in this form are depicted 
in Fig.iTjand Fig.[ 8 ] They may be more useful to experimental¬ 
ists than Fig.|5]and Fig.| 6 ] 


5 Remarks on experimental challenges 

Let us briefly motivate why, in the above presentation, we have 
(i) considered transport data for molecular junctions placed in 
solvent and (ii) referred to SERS as the technique to reveal the 
substantial dependence on the molecular charge of the vibra¬ 
tional properties predicted by our calculations. 

(i) Given the impossibility of achieving a substantial re¬ 
duction of the molecular species (significant change in the 
molecular charge) by varying the source-drain bias I 4 , an effi¬ 
cient molecular orbital gating appears to be irreplaceable. The 
molecular junctions studied in ref. exhibit the most sub¬ 
stantial orbital gating effect known to date in “dry” molecu¬ 
lar electronics. We checked by calculations similar to those 
presented above that, for those junctions, changes in molecu¬ 
lar charges do not exceed a few percent. This fact can actu¬ 
ally be understood intuitively by examining, e.g., Eig. S7 of 
the supplementary information of ref. 1^ which shows curves 
that do not exhibit any peak; a substantial reduction (oxida¬ 
tion) can only be obtained if the gate potential Vq utilized 
in experiments^ sample sufficiently extended portions of the 
I — Vg transfer characteristics comprising a maximum (which 
is the counterpart of the maximum of the I — 77 -curve shown 
in Fig.©. 

In principle, with a molecular orbital gating in “dry” en¬ 
vironment substantially improved well beyond the present 
achievements,— changes in molecular vibrational properties 
comparable to those presented above could also occur in nano¬ 
junctions placed in vacuum, as demonstrated by the data in 
Table S3. But, according to the state-of-the-art in the held, a(n 
almost) complete reduction of the molecular species (chang- 
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Fig. 7 Dependence on the current I of the changes in frequencies 
and Raman scattering activities (in A^/a.m.u.) in variable bias mode 
for several representative vibrational modes given in the legend. The 
V*-range sampled in experiment— is indicated by the two black 
vertical dashed lines. 

ing the molecular charge by an entire electron) can only be 
reached via the efficient electrolyte gating discussed in the 
present paper. 

(ii) It is hard to conceive that lETS (inelastic electron tun¬ 
neling spectroscopy) ^^"^4 can be utilized for the present pur¬ 
pose. Rather than SERS, lETS is the choice of experimen¬ 
talists for revealing vibrational effects in the charge trans¬ 
port through molecular junctions,However, it requires 
cryogenic temperatures. Obviously, such conditions are to¬ 
tally inappropriate for molecular junctions immersed in sol¬ 
vents. The fact that SERS can be applied at room tem- 
perature^ti^ii2iS is an important reason for advocating this 
method. 

The experimental demonstration of the SERS-related ef¬ 
fects under full electrochemical control discussed above is still 
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Fig. 8 Dependence on the current / of the changes in frequencies 
and Raman scattering activities (in A^/a.m.u.) in constant bias mode 
(variable rj) for several representative vibrational modes given in the 
legend. 


pending, which may be related to certain experimental chal¬ 
lenges, which we briefly mention below. 

An obvious challenge is the strength of the SERS sig¬ 
nal. Valuable insight on the SERS enhancement factor can 
be gained from 3D-EDTD (three-dimensional finite-difference 
time-domain) theoretical simulations, as shown in earlier 
work.2iL2i^ In the absence of any experimental SERS infor¬ 
mation on the EC-STM junctions envisaged, a ETDT simula¬ 
tion would be too speculative and will not be attempted be¬ 
low. In view of its well-known critical dependence on the 
incident light polarization, surface morphology, nanogap con¬ 
figuration and size, an estimate of the SERS intensities can be 
realistic only if a sufficiently detailed experimental character¬ 
ization is available. Still, we think that the following pieces 
of experimental work on related nanosystems can be taken as 
positive signals that a combined transport-SERS study in EC- 
STM junctions like those considered in this paper is feasible; 
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the extraordinary enhancement of the Raman scattering from 
pyridine by shell-isolated gold nanoparticles, ^ the application 
of in situ shell-isolated nanoparticle SERS (SHINERS)' * un¬ 
der electrochemical condition for an interfacial redox reaction 
using the same (viologen) molecule as presently considered, 
and the combined SERS and conductance measurements un¬ 
der the same experimental conditions for scanning tunneling 
microscope break junctions of 44BPY placed in aqueous so¬ 
lution.^ In the last case,- a Raman enhancement factor of 
3.8x10^ has been obtained. In view of the similarity between 
the presently considered system and that of ref. we do not 
expect a dramatic reduction of the SERS intensity preventing 
observability. 

Implementing SERS on single-molecule junctions^ repre¬ 
sented an important advance in molecular electronics. The 
so-called “fishing mode” STM is the key for Raman signal ac¬ 
cumulation. A shift of the Raman signal was observed and 
related to the stress applied to the molecule. Likewise, lETS 
(inelastic electron tunneling spectroscopy) studies on single¬ 
molecule junctions also indicated some shifts in vibrational 
modes under stress. So, stress-driven shifts could be a signif¬ 
icant noise factor for the rj- and 14-dependent Raman signals. 
However, while this may be an issue in general, and further 
joint experimental and theoretical studies are needed to clar¬ 
ify it, we believe that stress-driven shifts are less important 
for the specific case considered in this paper. As noted in the 
ESIf, the viologen core (44BPY) is not directly contacted to 
(gold) electrodes in the molecular junctions used in experi¬ 
ment, but rather via alkyl linkers. These linkers mitigate the 
electrodes’ impact on the 44BPY core. Therefore, we do not 
expect a substantial impact on the intra-core vibrational modes 
similar to that caused by the change in the molecular charge 
considered above. Related to this issue, it is noteworthy that 
even for the smaller 44BPY molecule {i.e., the molecule with¬ 
out alkyl(H-thiol) linkers) vibration frequencies computed with 
gold atoms attached at the two ends of the molecule and only 
at one end are almost equal; see Table S1 of ref. 

Obviously, neat curves like those presented in Fig. |6] and 
[8] represent a highly idealized description. In real measure¬ 
ments, one should rather expect blurry curves, similar to those 
of previous reports.^i2i2 In addition to the (albeit perhaps less 
pronounced) stress-driven shifts mentioned above, other dele¬ 
terious effects such as heatin g and photocurrents^“— 

are possible sources of noise. But just because such effects 
were also present in earlier successful SERS-transport stud¬ 
ies, there is no special reason to assume that they are so strong 
to prevent observability in the presently considered case. 

6 Conclusion 

To conclude, in this paper we have theoretically demonstrated 
that vibrational frequencies and Raman scattering intensities 


of a redox molecule embedded in an EC-STM single-molecule 
junction, like those already fabricated,— can be efficiently 
controlled especially by varying the overpotential (“gate” po¬ 
tential). As widely accepted, for SERS observability, the elec¬ 
tromagnetic enhancement is decisive although the chemical 
(charge transfer) enhancement— may also play a role.— The 
present results indicate that, in addition to these, the charge 
state may also be significant. The information that can be 
gathered by SERS goes beyond the chemical identity and valu¬ 
able structural information of the wired molecule and contacts; 
the charge state of a molecule can also be probed via SERS. 
Table S1 shows that differences in Raman intensities of the di- 
cationic and cationic species can be of several orders of mag¬ 
nitude. E.g., reduction yields the Raman activity of mode 23 
is enhanced by a factor of ~125 while for mode 4 it is di¬ 
minished by a factor ~640. Using other solvents, the Raman 
activity enhancement can even be much larger (almost 10"*^), 
as illustrated by the results of Table S2 for modes 10 and 14 
in benzene. This is an effect much stronger than previously 
reported in two-terminal setups without solvents.— 

Therefore, in spite of nontrivial experimental challenges 
mentioned above, we are confident that correlating transport 
data and vibrational information acquired by SERS from cur¬ 
rent carrying single-molecule junctions under electrochemical 
control will open new avenues of research that are not accessi¬ 
ble for nanotransport in two-terminal setup. Practical molec¬ 
ular electronics and fundamental science can be main benefi¬ 
ciaries. 
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